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Abstract 

The generalized Dirac equation of the second order, describing 
particles with spin 1/2 and two mass states, is analyzed. The projec- 
tion operators extracting states with definite energy and spin projec- 
tions are obtained. The first order generalized Dirac equation in the 
20-dimensional matrix form and the relativistically invariant bilinear 
form are derived. We obtain the canonical energy-momentum tensor 
and density of the electromagnetic current expressed through the 20- 
component wave function. Minimal and non-minimal electromagnetic 
interactions of fermions are considered, and the quantum-mechanical 
Hamiltonian is found. It is shown that there are only causal propaga- 
tions of waves in the approach considered. 



1 Introduction 

One of the modern problems is to explain the number of quark and lepton 
generations, and their mass spectrum. The standard model (SM) of elec- 
troweak interactions contains many free parameters: rrie, m^, rrir, m^, rrid, 
rric, ms, rrit, mb, four mixing angles and three coupling constants. One of the 
ways to reduce the number of free parameters is to explore the Grand Uni- 
fication Theories (GUTs) with different gauge groups. At the GUTs scale, 
we have only one coupling constant, and masses of leptons and quarks are 
generated below the GUTs scale by the spontaneous gauge symmetry break- 
ing with the help of the Higgs mechanism. As a result of the symmetry 
breaking, in the framework of the GUT, the generations of fermions appear. 
Nevertheless, the Higgs Lagrangians also contain free parameters. Recent 
observations of neutrino oscillations show that SM should be extended, and 
possibly a fourth generation exists. The deeper insight into the dynamics of 
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the mass generation mechanism is needed to understand physics beyond the 
SM. 

We pay attention here on the Barut work pp, (see also [2 J who suggested 
a mass formula for e-, /x-leptons based on the generalized Dirac equation of 
the second order (see [3], describing particles with two mass states (see 
[H], [H] for the case of bosonic fields with two mass states). This equation 
may be considered as an effective one for partly "dressed" fermions, i.e. some 
radiation corrections are taken into account. This scheme represents the non- 
perturbative approach to quantum electrodynamics. 

The goal of this paper is to formulate the mentioned second order equation 
in the form of the first order generalized Dirac equation (FOGDE), and 
investigate it. 

The paper is organized as follows. In Sec. 2, we define the mass spectrum 
of the generalized Dirac equation of the second order and find the projection 
operator extracting the solutions in the momentum space for free particles. 
The FOGDE is derived in Sec. 3. The relativistically invariant bilinear form, 
the Lagrangian formulation, and spin projection operators are given. In Sec. 
4, we obtain the canonical energy-momentum tensor and the electromagnetic 
current density from the Lagrangian. Sec. 5 is devoted to the introduction 
of minimal and non-minimal electromagnetic interactions of fermions in the 
approach discussed. The quantum-mechanical Hamiltonian is obtained. We 
make a conclusion in Sec. 6. The system of units h = c = 1 is chosen and 
notations as in |7] are used. 

2 Spin- 1/2 Field Equation of the Second Or- 
der 

Let us consider the second order (in derivatives) field equation describing 
spin-1/2 particles introduced in jS] (see also generalizations in [S]): 

(cti'ji.d,, + a2dl + Kj ip{x) = 0, (1) 

where c?j, = d/dxi, = {d/dxm,d/d{it)), ip{x) is a bispinor. Repeated indices 
imply a summation. The Dirac matrices 7^ obey the commutation relations 
Ifilu + lul^i = '^Sfj.u- As physical values depend only on two parameters, it is 
convenient to introduce new variables n/ai = m, 02/01 = —a/m, where m is 
a parameter with the dimension of the mass, and a is a massless parameter. 
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With these arrangements, Eq. (1) is rewritten as 

(7.5.-^a; + m)V'(^) = 0. (2) 
In momentum space Eq. (2) for positive energy is given by 

^ip + ^ + m^ = 0, (3) 

where p = 7^p^. One has to make the replacement p —p when considering 
antiparticles. Applying the operator {—ip + ap^/m + m) to both sides of Eq. 
(3), we obtain the equation 

aV + (2a + 1) + m^] V'(p) = 0, (4) 

where = — Po- order to have a non-zero wave function in Eq. (4), 
the diagonal elements need to be zero: 



2 2 I 2a + 1 ± V4a + 1 ^ 
P--m[ 1- (5) 



Eq. (5) defines the squared masses of fermions: 

m, = m ^ J, m2 = m ^ J. (6) 

From Eqs. (6), one finds the restriction on the parameter a: a > — 1/4. At 
a = —1/4 both masses are equal, mi = m2. We obtain from Eqs. (6) masses 
of spin- 1/2 particles 

m,^-m^ Ya J' = [ Ya J" 

In addition, we have a solution with opposite signs of masses. The signs of 
masses in Eqs. (7) are chosen to be positive for the negative parameter a. 
The negative mass solutions of Eqs. (6) can be incorporated in a description 
of antiparticles. When the parameter a approaches to zero, and the param- 
eter m is fixed, masses of the fermionic states, mi, m2, become unlimited. 
Therefore, the limit a = should be made in Eqs (2), (3), and one recovers 
the Dirac equation which describes fermions with the mass m. 
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Eqs. (7) lead to the mass formula mi + m2 = —m/a = axja^ obtained 
in 121 . The authors jS] interpreted the state of fermion with the mass mi as 
an electron, and the state of fermion with the mass m2 as a muon. 

Now we construct the projection operators extracting solutions of Eq. 
(3). It is easy to verify that projection operators (see the general method of 
projection operators in 0) 

ap^ + vn} =F imp 
^ " 2(ap2 + m2) ' 

which obey the necessary equation = A±, are solutions of equations 

I ±ip + — + m ) A± = 0. 
y my 

The projection operator A+ corresponds to positive energy, and A_ - to 
negative energy of particles. 



3 First Order Generalized Dirac Equation 

In order to formulate the FOGDE, we introduce the 20- dimensional function 

^{x) = {Mx)} = ( ) (M^) = -^OM^)), (9) 

where A = 0, fi. The function in Eq. (9) is the direct sum of a bispinor 
ip{x) and a vector-bispinor ip^{x). Therefore, the wave function realizes 
the reducible representation [(1/2,0) © (0, 1/2)] © {(1/2, 1/2) © [(1/2,0) © 
(0, 1/2)]} of the Lorentz group. 

Let us introduce the elements of the entire algebra e^'^ ^Ul with the 
properties 

[£ ) = OmaOn,b, e ' e ' = OabS , (10) 

y4, 5, M, = 0, 1, 2, 3, 4. The e^'^ represent 25 basis elements correspond- 
ing to the different values of the superscripts M, A^, and they are 5x5 matrices 
whose elements are labelled by the subscripts. The elements of the matrix 
^M,N consist of zeros and only one element is unity, where row M and column 
A^ cross. 
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With the help of Eqs. (9), (10), Eq. (2) can be rewritten in the form of 
the first order equation 

9, (e'^'O + a£°'" + V) ^s(x) +m + ^^(x) = 0, (11) 

where we imply that 7-matrices act on the bispinor subspace. All repeated 
indices such as in Eq. (11) imply a summation even in the superscripts of 
the elements of the algebra. After introducing 20-dimensional matrices 

a, = (e'^-o + a£°'") + ® 7,, (12) 



where unit four-dimensional matrix /4 acts on bispinor subspace and unit 
five-dimensional matrix 1 = = e^'^ + e^'^ acts on scalar-vector subspace, 
Eq. (11) takes the form of the FOGDE: 

[a^du + m) ^{x) = 0. (13) 

Eq. (13) is convenient for investigations of spin-1/2 fermions with two mass 
states. According to the general theory ^T], the masses of particles described 
by the first order relativistic wave equation (13) are m/Aj, where A, are the 
eigenvalues of the matrix a^. In our case the matrix 0:4 = (e^'*^ + ae^'^) ® 
I4 + ® 74 obeys the minimal matrix equation 

al-{l + 2a)al + = 0. (14) 

One can obtain from Eq. (14) four eigenvalues of the matrix a^. 

±A„ A.= -'-f^ . A.= -' + f^ . (15) 

It is easy to verify that masses of fermion states of 20-dimensional matrix 
Eq. (13): nii = m/Ai, 777,2 = '^/A2 coincide with those in Eq. (7). Two 
additional masses possess opposite signs. 

To build the spin operators, one needs the generators of the Lorentz group 
in the 20-dimensional representation. Such generators are given by ^21, QH] 

J,u = j!!J^h + h®J^^J'\ (16) 

J(l) = ^M,. _ (17) 
J^l^'^ = \il,lu-lul,). (18) 
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Operators J^^J , J^y^-* are the generators of the Lorentz group in four-dimensional 
vector and bispinor spaces, respectively. Generators (16)-(18) obey the com- 
mutation relations 

The antisymmetric parameters of the Lorentz group Umn {Tn,n = 1,2,3) are 
real, and are imaginary in our metric. The commutation relation 

["A, Jfiu] = hfiOiu - ^\yOl^ (20) 

holds and guarantees the relativistic form-invariance of Eq. (13). 

Now we consider operators of the spin projections on the direction of the 
momentum p: 

i 

= ^'^\ — \^abcVaJbc-, (21) 
^IPI 



where |p| = ypf + pi + vi- Using Eq. (20), it is not difficult to prove that 
the operator of Eq. (13) in the momentum space {ia^p^ + m) commutes 
with the spin operator (21): [ia^p^ + m, cTp] = 0. Therefore, the operators 
{idfiPii + ^) and dp have the common eigenfunction in the momentum space. 
It is not difficult to verify that the minimal matrix equation 



holds. With the aid of the method jUl, we find the projection operators 

1 / , U / 2 9 



P..,^ = =F5 ± ["i - jj (22) 
which extract spin projections ±1/2, so that 

crpP±i/2 = ±-P±i/2- 

The operators P±i/2 acting on the arbitrary 20-dimensional function \E'o(p) 
produce the eigenfunction '^±i/2{p) = -P±i/2^o(p) of the operator ap-. o"p \l'±i/2 ■ 

±(l/2)^±l/2. 

To obtain the Lagrangian, one has to find the Hermitianizing matrix rj in 
our 20-dimensional representation space which obeys the relations |lTj 

rjam = -o^mV: = (^tv ("^ = 1^ 2, 3). (23) 
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Such a matrix exists, and is given by 

V 



74- 



Using Eqs. (23), we obtain from Eq. (13) the equation 



^(x) (a^'S, 



m 



0, 



(24) 



(25) 



where ^(x) = '^'^{x)ri, and \l/"^(x) is the Hermitian-conjugate wave function. 
The derivative 9 ^ in Eq. (25) acts on the left-standing function. The 
relativistically invariant bihnear form is given by '^{x)'${x) = '^~^{x)ri'${x). 
Thus, the Lagrangian reads: 



1 r- 



£ = -- [^(x) [af^df, + m) ^(x) - ^(x) [a^^d^- m) ^(x) 



(26) 



Variation of the action S = J Cd^x with respect to the independent fields 
\Ef(x), ^E'(x), corresponding to Lagrangian (26), gives Euler-Lagrange equa- 
tions (13), (25). 



4 The Energy- Momentum Tensor and Elec- 
tromagnetic Current 

The energy-momentum tensor can be found with the help of the standard 
procedure, using the general formula [H] 

dC — dC 

= .T,r xN ^^^(^) + 5.^(a;)-7-^--Y " ^>^uC. (27) 
9 09;,^ X d(d^^ixy 



With the help of Eq. (27), one may obtain from the Lagrangian (26) the 
canonical energy-momentum tensor as follows 

T^u = \ (9.^(x)) a^^{x) - ]^{x)a^d,^{x). (28) 

We took here into consideration that £ = for functions \E'(x), \E'(x) obeying 
Eqs. (13), (25). Canonical energy-momentum tensor (28) is conserved but is 
not symmetric. It can be symmetrized ^ by adding the divergence term to 
the Lagrangian (26) which does not change Euler-Lagrange equations (13), 
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(25). With the help of Eqs. (9), (12), one may obtain from Eq. (28) the 
expression 

(df,d^^P{x)) 4j{x) + (d^ip{x)) df^tpix) - ^ip{x)d^d^ip{x). 
2m ^ ^ 2m ^ ' 2m 

The first two terms in this equation correspond to the ordinary Dirac equa- 
tion. The energy density and the momentum density are given by £^ = T44, 

The conserved electric current density may be obtained by the relation- 
ship Hi] 

^'<^-'^'(*'^->^kJm)-«*'4 '''' 

Substituting Eq. (26) into Eq. (29), we get the electric current density 

j^{x) = {^{x)a^^{x). (30) 

From Eqs. (13), (25), one may prove the conservation of the four- vector 
current density: d^j^{x) = 0. The canonical energy-momentum tensor (28) 
and electric current density (30) are expressed through 20-component wave 
function \I^(x). Using Eqs. (9), (12), we obtain from Eq. (30) the electric 
current density only in terms of the bispinor iplx): 

jf^i^) = -ii^{x)-f^ij{x) + —'il:{x)df,ij{x) - — {df,ip{x)) ip^x). 

It follows from this expression that the current (30) is related to Barut's 
original anomalous magnetism. So, the jfj.{x) includes the usual Dirac current 
as well as convective terms. 



5 Electromagnetic Interactions of Fermions 

It follows from Eq. (14) that the inverse matrix aj^ exists, as there are no 
zero eigenvalues of the matrix 0:4. One may find from Eq. (14) the inverse 
matrix 



2a + 1 

— 0:4 ■ 



a4 - -i- (31) 
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This indicates that all components of the wave function "^{x), Eq. (9), 
are canonical, and contain time derivatives. Thus, there are no subsidiary 
conditions in the first order wave equation (13), and it can be rewritten 
in the Hamiltonian form. The minimal interaction with electromagnetic 
field in the first order equation (13) can be obtained by the substitution 
dfj, ^ D^ — d^ — leAn [A^ is the four-vector potential of the electromagnetic 
field). We also introduce non- minimal interaction by adding terms linear 
in the strength of the electromagnetic field = dfj^A^, — dj^A^. So, we 
postulate the matrix equation 

i 

a^D^j, + - (ko-Po + i^iPi) OLixvT^v + m 
where the projection operators Pq; P\i are given by 

Po = ® ^4, ^1 = e^'^ ® h, (33) 

and 

= ® (7mO'. - 7.7m) + ® 7. - ® 7m (34) 

The projection operators Pq, Pi extract scalar, vector subspaces and obey the 
relations: Pq = Pq, Pf = Pi, Pq + Pi = 1. It is easy to verify that Eq. (32) is 
form-invariant under the Lorentz transformations. We have introduced two 
parameters Kq, Ki which characterize anomalous electromagnetic interactions 
of fermions. It should be noted that the relativistic invariance allows us to 
introduce only two parameters because there are two subspaces corresponding 
to the projection operators Pq and Pi. 

With the help of Eqs. (9), (10, (12), we obtain the tensor form of Eq. 
(32): 

{'y„D^ + iKo7M7i/-^Mi' + '^i^) + (o-^m + ii^olu^un) ^n{x) = 0, (35) 

{D/j, + ijj{x) + {m5fj_„ miaT^i,) V'i.(x) = 0. (36) 

It should be stressed that Eq. (36) is an analog of the gradient equation (see 
Eq. (9)) in the case of interacting fields. System of Eqs. (35), (36) defines 
bispinor il^{x) as well as vector-bispinor ipy{x). Eq. (36) may be considered as 
a matrix equation for vector-bispinor il^^ [x) . After finding the inverse matrix 



*(x) = 0, (32) 
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{mS^jj + miaJF^j,)"^, one can express the il^u{x) from Eq. (36) and replace 
into Eq. (35) to get an equation for bispinor ip{x) describing the interaction of 
fermions with electromagnetic fields. This equation includes the interaction 
with the anomalous magnetic moment of particles. Eqs. (35), (36) can also 
be applied for phenomenological descriptions of composite fermions. 

Let us find the quantum-mechanical Hamiltonian from Eq. (32). Eq. 
(32) can be rearranged as 



iaidt'^{x) 



ttaDa + m + eAoa4+ 



I 



^[x). (37) 



One can obtain from Eq. (37) the Hamiltonian form of the equation 



U 



% 

aaDa + m + eA^a^ + - (kq-Pq + i^iPi) 0!,,u^/,u 



(38) 



where the matrix aj^ is given by Eq. (31). 

To investigate the consistency problem, in accordance with the method 
[T3] . we replace the derivatives in Eq. (32) by the four- vector to get 
the characteristic surfaces. One can obtain the normals to the characteristic 
surfaces, n^, by solving the equation 

det (a^n^) = 0. (39) 

Using the frame of reference where ra^ = (0,0, 0,714), Eq. (39) becomes 
det(a4n4) = 0. It follows from Eq. (14) that there are no zero eigenvalues of 
the matrix 04. Therefore, equation det(a4ri4) = possesses only the trivial 
solution n4 = 0, that indicates: there are only causal propagations of waves 
as in the Dirac equation, and no superluminal speed of waves. 

6 Conclusion 

The projection operators extracting states with definite energy and spin pro- 
jections were constructed for the generalized Dirac equation of the second 
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order, describing particles with spin 1/2 and two mass states. We have ob- 
tained the FOGDE in the 20-dimensional matrix form, the relativistically 
invariant bihnear form, and the Lagrangian. The FOGDE is convenient for 
different apphcations. The canonical energy-momentum tensor and density 
of the electromagnetic current have expressed through the 20-component 
wave function. We have introduced two parameters characterizing non- 
minimal electromagnetic interactions of fermions. As a particular case, they 
include the interaction of the anomalous magnetic moment of particles. The 
quantum-mechanical Hamiltonian has obtained and the Hamiltonian form of 
the equation was given. It was shown that there are only causal propagations 
of waves in the approach considered. 

In order to get the third generation of fermions, one needs to introduce a 
third-order (in derivatives) generalized Dirac equation. Using the technique 
described, such equation can be reduced to first-order equation such as Eq. 
(13). There is not a natural hmit to these procedures, corresponding to a 
finite number of generations in the approach considered. 

The scheme investigated can be applied for a consideration of the two 
flavour generations of leptons and quarks. Possibly, Eqs. (35), (36) may be 
considered as effective equations for electromagnetic interactions of hadrons. 
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